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Abstract 

There are two type of scales present simultaneously in the space-like as well 
as in the time-like directions in a model-class describing a cylindrically sym- 
metric, finite, three-dimensionally expanding boson source. One type of the 
scales is related to the finite lifetime or geometrical size of the system, the 
other type is governed by the rate of change of the local momentum distri- 
bution in the considered temporal or spatial direction. The parameters of 
the Bose-Einstein correlation function may obey an Mf scaling, as observed 
in S -|- Pb and Pb + Pb reactions at CERN SPS. This M^-scaling may imply 
that the Bose-Einstein correlation functions view only a small part of a big 
and expanding system. The full sizes of the expanding system at the last in- 
teraction are shown to be measurable with the help the invariant momentum 
distribution of the emitted particles. A vanishing duration parameter can also 
be generated, with a specific Mt dependence, in the considered model-class. 
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I. INTRODUCTION 



The method of intensity interferometry has recently become a widely used tool for de- 
termining the space-time picture of high energy heavy ion collisions. Originally, the method 
was invented to measure angular diameters of distant stars. The objects under study 
were approximately static and the length-scales astronomical. In principle the same method 
is applied to measure space-time characteristics of high energy heavy ion collisions, where 
the objects are expanding systems, with life-times of a few fm/c (10~^^ sec) and length-scales 
of a few fm (10~^^ m). 

In the case of high energy heavy-ion collisions intensity interferometry is pursued to infer 
the equation of state and identify the possible formation of a transient Quark-Gluon Plasma 
state from a determination of the freeze-out hyper-surface, as scanned by the Bose-Einstein 
correlation function (BECF), see e. g. the contributions of the NA35, NA44 and WA80 
collaborations in ref. 0,0]. For introduction and review on Bose-Einstein correlations see 
refs. Non-trivial effects arising from correlations among space-time and momentum- 

space variables were studied in refs. 0. 

The recent '^^S' + ^^'^Pb reactions at 200 AGeV laboratory bombarding energy resulted 
in a non-expected, symmetrical BECF-s if measured in the LCMS, the longitudinally 
comoving system of the boson pairs The longitudinal radius parameter was shown to 
measure a length-scale, Rl oc Xj^JnTf, introduced in ref. for an infinite, longitudinally 
expanding Bjorken tube. The side radius parameter was thought to measure the geometrical 
radius and the out component to be sensitive to the duration of the particle freeze-out 
times [0,0. All radius component parameters turned out to be equal within the experimental 
errors. Although this might he just a coincidence, in this work we show that such a behavior, 
valid in a certain nit interval, may also be a natural consequence of a cylindrically symmetric 
three-dimensional hydrodynamic expansion. In this case the local temperature, the gradients 
of the temperature distribution and the flow-gradients generate 'thermal' length-scales in all 
these space-like directions. Changes in the local temperature during the particle emission 



2 



induce a temporal scale, the thermal duration. Recently it became clear that the parameters 
of the BECF-s measure the lengths of homogeneity P,p!0|-p!2| which in turn were shown to 



be expressible in terms of the geometrical and the thermal lengths ||13|,PJT2 . 

We shall derive here general relationships among the functional forms of the BECF-s as 
given in the laboratory (LAB) frame, the LCMS frame and the longitudinal saddle-point 
system (LSPS) in which the functional form of the BECF-s turns out to be the simplest one. 

A new class of analytically solvable models is introduced thereafter, describing a three- 
dimensionally expanding, cylindrically symmetric system for which the geometrical sizes and 
the duration of the particle emission are finite. In this class of the models there are two 
length-scales present in all directions, including the temporal one. The BECF is found to be 
dominated by the shorter, while the momentum distribution by the longer of these scales. 
The interplay between the finite "geometrical scales" of the boson-emitting source and the 
finite "thermal scales" shall be considered in detail. 

II. FORMALISM 

Both the momentum spectra and the BECF-s are prescribed in the applied Wigner- 
function formalism |TB|,|^. In this formalism the BECF is calculated from the two-body 



Wigner-function assuming chaotic particle emission. In the final expression the time- 
derivative of the (non-relativistic) Wigner function is approximated |l7|JT5[| by a classical 
emission function S{x;p), which is the probability that a boson is produced at a given 
X = (t,r) = (t, ) point in space-time with the four-momentum p = {E,p) = 

{E,px,Py,Pz). The emission function has been related to the covariant Wigner-transform of 



the density matrix of pion sources in refs. [|T^,§] and most recently in ref. |T^, where the 



relation of Wigner-function formalism to the covariant current formalism has also been 
clarified. The (off-shell) two-particle Wigner functions shall be approximated by the off-shell 
continuation of the on-shell Wigner-f unctions |]r5| , p!3| ,p| . The particle is on the mass shell. 



= — p^. Please note the difference between x indicating a four- vector in space-time 
and the script-size x which indexes a direction in coordinate space. 



A useful auxiliary function is the Fourier-transformed emission function 

S{Ak; K)= j (fx S{x] K) exp{iAk ■ x), (1) 

where 

P1+P2 

Ak = p^-P2, K= (2) 

2 

and Ak ■ x stands for the inner-product of the four-vectors. Then the one-particle inclusive 
invariant momentum distribution (IMD) of the emitted particles, A^i(p) is given by 

E da 

N,{p) = S{Ak = 0;K = p) = , (3) 

cTtot dp 

where atot is the total inelastic cross- sect ion. This IMD is normalized to the mean multi- 
phcity (n) as 

/ -iVi(p) = (n). (4) 
E 

In the present paper effects arising from the final state Coulomb and Yukawa interactions 
shall be neglected. The two-particle BECF can be calculated from the emission function 
with the help of the well-established approximation 

iV2(pi,p2) \S{Ak;K)\' 
C{Ak;K) = ^1 + , (5) 

(n(n-l))iVi(pi)iVi(p2) I S{0;K) \' 



utilized also in ref . , see ref . |jT5| for further details. The corrections to this expression are 



known to be small [|T^. Note that among the eight components of Ak and K only six are 
independent due to the two constraints Pi = P2 = These constraints can be formulated 
alternatively as Ak ■ K = and = m? — Thus the two-particle BECF depends 

on the off-shell emission function, which we approximate by the off-shell continuation of the 
on-shell emission functions. 

A similar but not identical approximation used by several authors is to replace S{Ak] K) 
by S{Ak; K') where off-shell K is changed to on-shell K' . The latter mean momentum is 



defined to be on-shell as K'^ = w? — K'^ where K' = K = (pi + P2)/2. The differences 
between these two approximation schemes are of (9(Afc^/m^). The above two approximation 
schemes coincide in the AA;^ hmit where the Bose-Einstein correlations are maximal 
Since we shall make use of the ISk-K = ^ constraint which is exact only if the K four- vector 
is off-shell, we shall approximate the off-shell emission function in eq. ( |^) with the off-shell 
continuation of the on-shell emission function. 

III. GENERAL CONSIDERATIONS 

We model the emission function in terms of the longitudinally boost-invariant vari- 
ables: the (longitudinal) proper-time, r = yr^ — rf, the space-time rapidity r] = 
0.5 In ( (t + Vz) I {t — Vz)), the transverse mass nit = \J E'^ ^ vl ^ind the momentum-space 
rapidity y = 0.5 In ( (i? + p^) / — p^) ). In the transverse direction, the transverse radius, 
rt = \Jr1 + Ty is introduced. We have 

t = rcosh(?7), = r sinh(?7). (6) 

For systems undergoing a boost-invariant longitudinal expansion, the emission function may 
be a function of boost-invariant variables only. These are r, rx,ry, Px,Py and rj — y. How- 
ever, for finite systems the exact longitudinal boost-invariance cannot be achieved and the 
emission function becomes a function of rj — yo too, where yo stands for the mid-rapidity. 
Approximate boost-invariance is recovered in the mid-rapidity region only, where terms pro- 
portional to rj — yo can be neglected. Thus for finite systems undergoing a boost-invariant 
longitudinal expansion the emission function can be given in terms of these variables as 

S{x; K) d^x = S^{t, 1], fx, Ty) dr Todr] dr^ dry. (7) 

Here we introduced the constant tq in front of drj due to dimensional reasons and included the 
Jacobian from the d'^x to the dr drj dr^ dry variables into the emission function S^,{t, t], r^, ry) . 
The subscript * indicates that the functional form of the emission function is changed with 
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the change of the variables. Further, dependences on the mean momentum K as well as 
on the mid-rapidity yo are also indicated with the subscript The effective, momentum- 
dependent parameters of the emission function S^,{T,r],rx,ry) shall also be indexed with ^, 
in the forthcoming. The subscript s stands for the point where the emission function is 
maximal (we assume that S{x; K) has only one maximum for any values of K). We do not 
assume at this point whether the function — In S{x\ K) is expandable into a (multi-variate) 
Taylor series fl^ around its unique minimum at the saddle point Xs or not, merely we assume 
that the Fourier-transformed S'(Afc; K) exists. See the Appendix for a clarifying example. 
We suppose, however, that the Fourier-transformed S{Ak; K) can be evaluated in terms of 
the r and rj variables in the small Ak region relevant for the analysis of the BECF-s. This 
is possible if the region around Xs{K), where the Fourier integrals pick up the dominant 
contribution from, is sufficiently small so that within this region the r and rj dependence of 
t and Tz can be linearized as 

t ~ r cosh.[r]s] + {v - Vs)'rs smh.[r]s], (8) 
~ r sinh[?75] + {i] - r]s)Ts cosh[r7^], (9) 

with negligible second-order corrections. This condition is fulfilled if the characteristic sizes 
Ar* and At] of the considered region around Xs{K) satisfied Ar^ << r| and Ar]l << 1. 

The principal directions for the decomposition of the relative momentum at a given 
value of the mean four- momentum K are given as ||T8|,[7[|: the out direction is parallel to 



the component of K, which is perpendicular to the beam, indexed with out, the longitudinal 
or long direction is parallel to the beam- axis r^, this component of the relative momentum 
is indexed with l, and the remaining direction orthogonal to both longitudinal and out is 
called the side direction, indexed with side- Thus the mean and the relative momenta are 
decomposed as = {Kq, Kout, 0, Kl) and Ak = {Qo, Qout, Qside, Ql)- 
Since the particles are on mass-shell, we have 

Q = K-Ak = KoQo - KoutQout - KlQl. (10) 
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Thus the energy difference Qo can be expressed as 

Qo = PoutQout + PlQl, (11) 

where we have introduced the longitudinal and the outward component of the velocity of 
the pair, Pl = Kl/Kq and Pout = Kout/Ko, respectively. These relations become further 
simplified in the LCMS, the longitudinally co-moving system, introduced first in ref. |J7|. The 
LCMS is the frame where Kl = thus /3l = 0. We also have Pout = Pt where t stands for 
transverse, e.g. = + and rrit = + p1+ pI- Note that the relation Pout = Pt is 
independent of the longitudinal boosts, but both sides of this equation transform hke I/Kq. 

Let us express the Fourier integrals in terms of the r and r] variables in the laboratory 
reference frame (LAB), utilizing eqs. (HJ^). The results in LCMS can be obtained from the 
more complicated results in the LAB frame by the substitution Pl = and Pout = Pt- To 
simplify the notation, let us rewrite 

Ak-X = Qot - QoutTx - QsideTy - Q lT z - QrT - QoutTx " QsideTy - Qr,rs{r] - 1]^), (12) 

utilizing the linearized eqs. (|§,|^). We have introduced the coefficients of the r and the 
Ts{r] — i]s) as new variables given by 

Qt = Qo cosh.[i]s] - Ql sinh[?73] = {PtQout + PlQl) cosh\ris] - Ql sinh[?7^], (13) 
Qr, = Ql cosh[r)s] - Qo sinhf?]^] = Ql cosh[?7s] - {PtQout + PlQl) sinhf?]^]. (14) 

From these relations it follows that 

I S{/S.k',K) I I Sif^QriQviQout^Qside) \ 

C{Ak;K)c^l + -1 + . (15) 

\S{0;K)\^ I 5,(0,0,0,0) |2 

Note that this expression contains a four- dimensional Fourier-transformed function, and 
among the four variables Q^Qr), Qout and Qside only three are independent due to eq. (^Ij). 
Note also that at this point the BECF may have a non-Gaussian structure, and its depen- 
dence on its variables does not factorize. The main limitation of the last approximation in 
eq. (^) is that it is valid only for systems with small lengths of homogeneity, At^ << 
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and JS.r]1 << 1. As we shall see in the forthcoming, this gives a lower limit in mt for the 
applicability of the simple analytic results for a certain class of emission functions. 

IV. THE CORE/HALO MODEL 

If the system under consideration consists of a core characterized by a hydrodynamic 
expansion and small regions of homogeneity, and a surrounding halo of long-lived resonances, 
then the above general expression can be further evaluated if the halo is characterized by 
sufficiently large regions of homogeneity. Indeed, the long lived resonances may decay in a 
large volume proportional to their lifetime, and the decay products are emitted with a given 
momentum distribution from the whole volume of the decay. 

The key point is the following: Let us consider an ensemble of long-lived resonances 
with similar momentum, emitted from a given small volume of the core. The momentum 
distribution of the decay products of these resonances will be similar to each other, inde- 
pendently of the approximate position of the decay. Now the approximate position of the 
decay is randomly distributed along the line of the resonance propagation with the weight 
P{t) oc ex.^[—mres^rest/Eres)- Thus the dccay products will be emitted with the same 
momentum distribution from a volume which is elongated along the line of resonance prop- 
agation, given by Vd^cay — I Pres I I ij^res^ res) 1 whcrc Aq is the initial transverse size of 
the surface through which the resonances are emitted with a momentum p^es approximately 
at the time of the decay of the core, Tg. 

Thus the halo of long-lived resonances is characterized by large regions of homogeneity. 
(In case of the pionic halo the dominant long-lived resonances are u!,r],rj' and K^, all with 
life-times l/F^es greater than 20 fm/c). If the emission function is a sum of the emission 
function of the core and the halo, 

S(x- K) = 5,,,(t, 77, r^, ry) + Sh{x- K) , (16) 

and the Fourier-transformed emission function of the halo is sufficiently narrow to vanish 
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at the finite resolution of tlie relative momentum Ak in a given experiment, then one can 



show El that 



iVi(p) =iVi,,(p)+iVi,,(p), (17) 

I q(^Qt-j Qri) Qouti Qside) \ 

C{/\k]K) = l + K , (18) 

I 5,,e(0, 0,0,0) |2 

where A'^i^j(p) indicates the number of particles emitted from the halo or from the core for 
i = h,c and the effective intercept parameter, 

2 



A* = A*(ir ~ p) 



(19) 



iVi(p) 

is the square of the ratio of the number of particles emitted from the core to the number 
of all the emitted particles with a given momentum p. This effective intercept parameter 
arises due to the finite relative momentum resolution and the comparably large region of 
homogeneity characterizing the halo part of the system. 

We would like to raise a warning flag here: The volume, which the decay products of the 
long-lived resonances of a given momentum are emitted from, is large only if the decaying 
resonances have | p | /{nires^res) » 1 fni/c. This in turn implies that the above simple 
picture may need further corrections for very low pt pions at rapidity y = 0. 

There is a gap in the life-time distribution of abundant hadronic resonances: l/Fp ^ 1.3 
fm/c, I/Ftv* — 0.56 fm/c, I/Fa — 1.6 fm/c and l/T^* — 3.9 fm/c, which life-times are 
of the same order of magnitude as the time-scales for re-scattering at the time of the last 
hadronic interactions. These lifetimes are also all a factor of 5 - 10 shorter than the life-time 
of the u meson, which is the long-lived resonance with shortest lifetime. Thus the decay- 
product of the short-lived resonances will mainly contribute to the core, which is resolvable by 
BEC measurements, while the decay-products of long-lived hadronic resonances will mainly 
belong to the halo, re-defined alternatively as the part of the emission function which is not 
resolvable in a given Bose- Einstein measurement. 
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V. CLASSES OF SIMPLE CORE FUNCTIONS 



/f the emission function of the core can be factorized, 

S*,cir, r], r^., Vy) = H^{t) G^{r]) /^r^., r^), (20) 

where H^{t) stands for the effective emission function in proper-time, G^{ri) stands for 
the effective emission function in space-time rapidity, and h{rx,ry) stands for the effective 
emission function in the transverse directions, then the expression for the BECF can be 
further simphfied as 

I -f^*(Qr) I I Gjfi^Qri) I I I*(^Q outi Q side) \ 

C{Ak;K) = l + X, . (21) 

iif,(o)na(o)n/*(o,o)p 

If the /*(r^., Ty) function is symmetric for rotations in the (r^., Vy) plane around its maximum 



point Vx^s then one may introduce Qt = y Qiide + Qout to find 

I H^{Qt-) I G^{Q.n) 1^ I h{Qt) r 
C{M-K) = l + K . (22) 

1^,(0) na(o) n/.(o) |2 

Such factorization around the saddle-point happens e.g. for the new class of analytically 
solvable models if certain conditions are satisfied, as discussed in the subsequent part. From 
the above expression it is clear that for this type of models the dependence of the BECF 
on the components of the relative momentum can be diagonalized with appropriate choice 
of the three independent components of the relative momentum. Note that the assumed 
existence of the Fourier-transformed distribution functions is a weaker condition than the 
assumption of the analytic form of the Fourier-transformed function, see Appendix for an 



example. Another example was given e.g. in ref. |T9[ for a H{t) distribution for which H{Q^) 
is not analytic function at = and | H{Qt) p does not start with a quadratic term. In 
mathematical statistics it is well known that the Fourier-transformed stable distributions 
are not analytic at Q = |2^. On the other hand, there are many physically interesting 
Gaussian models which correspond to the multivariate second order Taylor expansion of 
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the above general results, i.e. the analytic form of the corresponding Fourier-transformed 



function. The out-longitudinal cross term |12| has been recently discovered also in this 



context. To study the properties of the BECF let us apply a Gaussian approximation to the 
effective distribution functions as 

H^ir) oc exp(-(r - r,)V(2Ar2) ), (23) 
GM « exp(-(r/ - Vs f/{2Ar]l ), (24) 
hir^^Ty) oc exp(-( (r,. - r,,,)^ + {ry - ry^sf)/{2Rl) ). (25) 

Apart from the momentum- dependent parameters Ar^,, A?7* and the mean emission point 
may also be momentum- dependent in the above expression, = Ts{K)^ rjg = ris{K), r^^s = 
^x,s{K) and Ty^g = ry s{K). For the sake of simplicity we do not specify the normalization 
constants in eq. (^) since they cancel from the BECF which is given by 

CiAk; K) = l + X, exp{-QlAr^ - Qy.Aiil - QIRD- (26) 

This is a diagonal form of BECF-s for which the factorization property, eq. (pOf ) and the 
Gaussian approximation for the core, eqs. (p3|-p^) are simultaneously satisfied. In the present 



form of the BECF, there are no cross-terms among the chosen variables. Now, let us rewrite 



this form using the standard HBT coordinate system ||18[ to find 

C{Ak- K) = l + K exp(-i?LeQLe - RlutQlut - RlQl - 2Rl,^LQoutQL), (27) 

R-lide = R-ly (28) 

R-lut = Rl + ^Rluty (29) 

SRlut = f3K^osh'[vs]AT^ + smh'[v,]T^Ar^% (30) 

Rl = sinh[r/,] - coshir^.^VfAr/^ + (/^^ cosh[r/,] - sinhfr^.D^Ar^, (31) 
Rlut,L = (Acosh[r/,](/5LCOsh[r7,] - sinh[r/,]))Ar2 + 

(A sinh[77,] smh[r/,] - cosh[r7,]))r2 Ar/^. (32) 

This result is non-perturbative in terms of the variable rj^ and is valid in any frame. The 

main limitations of this result are the assumed Gaussian model-class, c.f. eqs. (p3| - p5D and 



11 



the assumed smallness of the emission region around Xs{K) so that t and dependencies 
could be hnearized in terms of r and rj. 

The above equations simphfy a lot in the LCMS system, where /5l = 0: 

SRLt = f3Kcosh%,]AT: + sinh'[vs]T^Avl), (33) 
Rl = cosh^ [r]s]r^Ar]l + sinh^ [r/,] Ar^, (34) 
Rlut,L = -Pt sinh[r/,] cosh[r/,](Ar,2 + r^Ar^l). (35) 

The life-time information Ar^ and the invariant measure of the longitudinal size along the 
Ts = const hyperbola, r^Arjl appear in a mixed form in the B^^,R\ and the Rl^t l source 
parameters even in the LCMS frame. The amount of these mixings is controlled by the 
value of Tj^^^'^^ . This relationship clarifies the physical significance of the rj^'^'^'^^ , the space- 
time rapidity of the maximum of the emission function in the LCMS frame: rj^^'^^^ is the 
cross-term generating hyperbolic mixing angle for cylindrically symmetric, finite systems 
undergoing longitudinal expansion and satisfying the factorization property eq. (pOD. If 
^LCMS = 0, no mixing of temporal and longitudinal components appear in LCMS. In some 
limited sense one may call rjs the cross-term generating hyperbolic mixing angle in any frame, 
because if = in a certain frame than cross-terms can be diagonalized away as follows. 

Let us define the LSPS, the longitudinal saddle point system, to be the frame where 
rjs = 0. Since rjs is a function of in a fixed frame, rjs = t]s{K), the LSPS frame may 
depend on K (e.g. on transverse mass of the pair). In the LSPS frame the out-long cross- 
term and the mixing of the temporal and time-like informations can be diagonalized. We 
have in LSPS 

SRlut = Pt^rl (36) 
Rl = T^Ar^l + PlATl (37) 
Rlut,L = PtPiArl (38) 



as follows from eqs. (|30| - p2D . Introducing the new variables Qq = PtQout + (^lQl and Qt 



Qout + Qlide '^^ obtain for the correlation function 
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C(Afc; i^') = 1 + A. expi-Ar^Q', - r^AT^tQi - RiQl). (39) 



From this relationship we also see that Qo{LSPS) = Qr,QL{,LSPS) = Qn, c.f. eq. (|26D . 

Let us study an expansion in terms of e =| Y — yo \ /At], where Y is the rapidity 
belonging to K the mean momentum of the pair, and Arj is the geometrical size of the 
expanding system in the space-time rapidity variable, satisfying Arj > Ar/*. It is obvious 
that in the LAB frame r]^^^ = Y + 0{e), since in the e — > limit we recover boost-invariance 
and the particle emission must be centered around the only scale: the rapidity of the pair. 
Similarly we see that r^fc-*^-^ = + C(e). It follows that the cross-term and the crossing of 
temporal and longitudinal information in the LAB frame is a leading order effect, 

SRLt = A'(cosh2[F]Ar,2 + smh'[Y]T^Avl) + 0(e), (40) 

Rl = + 0(e), (41) 

cosh^[y] 

sinh[y] 

RLt,L = -A r^^vl + 0{e). (42) 

cosh[F] 

On the other hand, the mixing of the temporal and longitudinal information is only next-to 
leading order in the LCMS according to eq. (§^), i.e. R^^^i^{LCMS) = + 0(e). However, 
if the I Y — Uq \« Ai] condition is not satisfied, the out-long cross-term might be large even 
in LCMS, as has been demonstrated numerically in ref. |^ . 



The cross-term generating mixing angle r^^ vanishes exactly in the LSPS frame, becomes 
a small parameter in the LCMS if | F — I /^V << 1 and becomes leading order in any 



frame significantly different from LSPS or LCMS. Thus we confirm the recent finding p2 
that the out-longitudinal cross-term can be diagonalized away if one finds the (transverse 
mass dependent) longitudinal rest frame of the source. 

Note, that cylindrical symmetry around the center of the particle emission as as- 
sumed by eqs. (|23| - |25D is a stronger requirement than the cylindrical symmetry of the 
emission function around the beam axis. This latter symmetry implies only that 
both the requirements S{t,r^,ry,r^; Kq.KouuKl) = S(t, -r^,ry,r^; Kq, -KouuKl) and 
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S{t ; Ko,Kout,KL) = S{t,r^, -ry,r^; Kq^Koui^Kl) should be simultaneously ful- 

filled. Thus cylindrical symmetry around the beam axis is compatible with a different 
Gaussian radius in the side and the out direction, 



I*{rx,ry) oc exp 



(43) 



with R^^x 7^ Cylindrical symmetry around the beam axis implies only that 

rr,^s{Ko,Kout,KL) = -rr,^s{Ko, -Koui.Kl) and ry^s{K) = 0. In the low transverse mo- 
mentum limit, when Kout = 0, the relations r^^siKout = 0) = and R^,^x{,Kout = 0) = 
R^:y{Kout = 0) also follow from cylindrical symmetry around the beam axis. If i?*^^. 7^ R^^^y 
at a given non- vanishing value of the mean transverse momentum Kout-, the generalized 
version of eqs. (^39) for the BECF reads as 



C(Afc; K) = l + K eM-^rlQl - ^^ItIQI - R%Ql^, - i?^,,gLe)- (44) 
In such a case, eqs. (|^,^) are also modified as 

^side — (45) 

-RL* = ^*,x + ^^ouf (46) 

This implies that the difference between the out and side radius parameters is not restricted 
by cylindrical symmetry around the beam axis to positive values only, since -R^Mt ~ ^side — 
5R?g^^ -f R^ x — Rl y which can also be negative if is sufficiently large However, 
cylindrical symmetry does imply that Rout = Rside in the K^ut limit. 

Up to this point, we have reviewed the properties of BECF-s without reference to any 
particular model, for some more and more limited classes of simple emission functions. We 
have obtained certain model-independent relations c.f. eqs. ( p!5| , |l8| , p2D which are valid for 
some non-Gaussian as well a Gaussian source functions. We have studied the relations 
between source parameters with a method which is non-perturbative in terms of r/^, but 
perturbative in terms of Ar/^ and Ar^/r^^. 
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Let us study the properties of an analytically solvable model-class in the subsequent 
parts. 



VI. A NEW CLASS OF ANALYTICALLY SOLVABLE MODELS 



For central heavy ion collisions at high energies the beam or axis becomes a symmetry 
axis. Since the initial state of the reaction is axially symmetric and the equations of motion 
do not break this pattern, the final state must be axially symmetric too. However, in order 
to generate the thermal length-scales in the transverse directions, the flow-field must be 
either three-dimensional, or the temperature distribution must have significant gradients in 
the transverse directions. Furthermore, the local temperature may either increase during 
the duration of the particle emission because of the re-heating of the system caused by 
the hadronization and/or intensive re-scattering processes or decrease because of the 
expansion and the emission of the most energetic particles from the interaction region. An 
example for such a time-dependent temperature was given e.g. by the solid line of Fig. 1. 
in ref . |2j 

We study the following model emission function for high energy heavy ion reactions: 



S{x]K) d^x 
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rrit cosh.{ri — y) exp 



K ■ u{x) fi{x 
+ ■ 



H{T)dT T^drj dr^dry. (47) 



V T{x) T{x) j 

Here g is the degeneracy factor, the pre-factor rrtt cosh(r7 — y) corresponds to the flux of the 
particles through a r = const hypersurface according to the Cooper-Fry formula and the 
four- velocity u{x) is 
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U[X] 
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(48) 



which describes a scaling longitudinal expansion with a linear transverse flow profile. The 
transverse flow is assumed to be non-relativistic in the region where there is a significant 
contribution to particle production. The local temperature distribution T{x) at the last 
interaction points is assumed to have the form 
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(49) 



T{x) To 

and the local rest density distribution is controlled by the chemical potential n{x) for which 



we have the ansatz 



= (50) 

T{x) To 2Rl 2Ari^ 

The parameters Rq and At] control the density distribution with finite geometrical sizes. 
The proper-time distribution of the last interaction points is assumed to have the following 
simple form: 

H{t) = ^ exp(-(r-ro)V(2Ar2))). (51) 

(27rAr2)(V2) 

The parameter At stands for the width of the freeze-out hypersurface distribution, i.e. the 
emission is from a layer of hypersurfaces which tends to an infinitely narrow hypersurface 
in the At — > limit. 

This emission function corresponds to a Boltzmann approximation to the local momen- 
tum distribution of a longitudinally expanding finite system which expands into the trans- 
verse directions with a transverse fiow which is non-relativistic at the saddle-point. The 
transverse gradients of the local temperature at the last interaction points are controlled by 
the parameter a. The strength of the flow is controlled by the parameter b. The parameter 
c = 1 is reserved to denote the speed of light, and the parameter d controls the strength of 
the change of the local temperature during the course of particle emission. 

Note that other shapes of the temperature profile lead to the same result if 1/T{x) starts 
with the same second order Taylor expansion around = Vy = 0. The physical significance 
of the transverse temperature profile is that it concentrates the emission of the particles with 
high transverse mass to a region which is centered around r^; = Tj, = and which narrows 
as the transverse mass increases. The Gaussian approximation to the inverse temperature 
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profile is thus a technical simplification only, other decreasing temperature profiles have 
similar effects as follows from the above picture. Similarly, the significance of the temporal 
changes of the temperature is that it creates different effective emission times for particles 
with different transverse mass, the Gaussian approximation is technical simplification. 

For the case o{a = b = d = Owe recover the case of longitudinally expanding finite 
systems as presented in ref. [^. The finite geometrical and temporal length-scales are rep- 
resented by the transverse geometrical size Rg, the geometrical width of the space-time 
rapidity distribution At] and the mean duration of the particle emission At. Effects arising 



from the finite longitudinal size were calculated analytically first in ref. p5[ in certain lim- 
ited regions of the phase-space. We assume here that the finite geometrical and temporal 
scales as well as the transverse radius and proper-time dependence of the inverse of the local 
temperature can be represented by the mean and the variance of the respective variables 
i.e. we apply a Gaussian approximation, corresponding to the forms listed above, in order 
to get analytically tractable results. We have first proposed the a = 0, 6 = 1 and d = 
version of the present model, and elaborated also the a = b = d = model corresponding 
to longitudinally expanding finite systems with a constant freeze-out temperature and no 
transverse fiow. Soon the parameter b has been introduced |T2| and it has been realized 



that the maximum of the emission function for a given mean momentum K has to be close 
to the beam axis, fulfilling ^ << tq, in order to get a transverse mass scaling law for the 



parameters of the Bose- Einstein correlation functions in certain limiting cases |26]. In this 



region around the beam axis, however, the transverse fiow is non-relativistic [|T2| even for 
the case 6 = 1 if this region is sufficiently small. Yu. Sinyukov and collaborators classified 
the various cases of ultra- relativistic transverse fiows [|lT|, [El], and introduced a parameter 



which controls the transverse temperature profile, corresponding to the a ^ b = case. We 



have studied |^ the model-class a^O, b^O,d = which we extend here to the d ^ 
case too. 

The integrals of the emission function are evaluated using the saddle-point method 
||8| JT0|JT^ . The saddle-point coincides with the maximum of the emission function, parame- 
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terized by (t^, 77^, r^..^^, r^^s)- These coordinate values solve simultaneously the equations 

dS dS dS dS 

— = — = = = 0. (52) 

dr drj dr^ dry 

These saddle-point equations are solved in the LCMS, the longitudinally comoving system, 
for Tj^^^^^ « 1 and r^^s « To- The approximations are self-consistent if | Y — uq \« 
l+/S.rfmt/TQ—/\rf and /3t << T^Arjl / (bRl) which for the considered model can be simplified 
as = Pt/^t << (o^ + b'^)/b/ max(l,a, 6). The transverse flow is non-relativistic at the 
saddle-point if Pt « (a^ + 6^)/6^/ max(l, a, h). We assume that Ar < tq so that the Fourier- 
integrals involving H{t) in the < r < 00 domain can be extended to the —00 < r < cxd 
domain. The radius parameters are evaluated here to the leading order in ,,/ tq. Thus terms 
of 0{rx,s/TQ) are neglected, however we keep all the higher-order correction terms arising 
from the non- vanishing value of rjg in the LCMS. We calculate both the radius parameters 
and the invariant momentum distribution in Gaussian saddle-point approximation. We shall 
discuss the hmitations of the saddle-point method after presenting these results on BECF 
and IMD. 



For the model of eq. ( |47|) the saddle point approximation for the integrals leads to an 
effective emission function which can be factorized similarly to eq. (pO|). Thus the radius 
parameters of the model are expressible in terms of the homogeneity lengths Ar^*,/?*, Ar* 
and the position of the saddle point rjs i.e. the cross-term generating hyperbolic mixing angle. 
The saddle-point in LCMS is given by = Tq, ryf'C'*^-^ = (y^ _ + Ar]'^{l/Ar]l - 1)), 

fx,s = PtbRl/ {tqAt]^) and Vy^^ = 0. Note that the space-time rapidity of the saddle-point 
^LCMS (jgpgnds on the boost-invariant difference yo — Y which can be evaluated in any 



frame. The radius parameters or lengths of homogeneity |§|,|T2[ are given in the LCMS by 



eqs. ( Kn - |29|j33| - |35|) , and we obtain 
1 1 1 



— = — + — cosh[r]^^^'% (53) 



T-)2 p2 p2 
^* 
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Ill 1 

= + cosh[r7f'^^^] , (54) 

/\r^l Ar]' Ar^l cosh^ [r/f ^^^^] 

111 

= + cosh[r;f ^^^] . (55) 



Arl Ar^ Ar| 



where the thermal length-scales are given by 

= -, (56) 

+ Mt 

T 

Aril = — , (57) 
Mt 

Ar| = . (58) 

d^Mt 

Here Mt = \J Kq — Kl is the transverse mass belonging to the mean momentum K. In the 

region of the Bose-Einstein enhancement, where the relative momentum of the pair is small, 
1 

Mt satisfies Mt = -{rrit^i + mt,2){l + 0{yi - 2/2) + 0{{mt^i - mt^2)/{m,i + m,2)))- Please 
2 

note the distinction between the subscripts for transverse direction, indicated by t, and the 
subscripts for the 'thermal' scales indicated by t- It is timely to emphasize at this point 
that the parameters of the Bose-Einstein correlation function coincide with the (rapidity and 
transverse mass dependent) lengths of homogeneity in the source, which physically can 
be identified with that region in coordinate space where particles with a given momentum 
are emitted from. The above relations indicate that these lengths of homogeneity for simple 
thermal models can be basically obtained from two type of scales in the framework of the 
saddle-point method. These scales have different momentum - dependence and are referred 
to as 'thermal' and 'geometrical' scales. 

In contrast to the homogeneity lengths which can be defined even without thermaliza- 
tion, the 'thermal scales' cannot be introduced without at least approximate local ther- 
malization. Thus the thermal scales originate from the factor exp(— p ■ u{x) /T{x)). They 
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measure that region in space-time where thermal smearing can compensate the change of 
the local momentum distribution which in turn is induced by either the gradients of the flow- 
field or the gradients of the temperature field. This is to be contrasted to the 'geometrical' 
scales, which originate from the exp(/x(x)/T(a;)) factor which controls the density distribu- 
tion. The geometrical scales can be interpreted as the regions in space-time where there 
is significant density to have particle emission. Obviously for locally thermalized systems 
both the geometrical and the thermal scales influence the regions of homogeneity and the 
smaller scale will be the dominant one. Since the four-momentum p is explicit in the factor 
exp{—p-u{x)/T{x)), and enters the 'geometrical' scales only through the momentum depen- 
dence of the saddle-point, the momentum dependencies for the 'thermal' and 'geometrical' 
scales shall be in general different from each other. Note also that in the above expression 
for A?7^ a third type of scale is also present in the term —1/ cosh^[?7s], which stems from the 
mtcosh[?7 — y] Cooper- Fry pre-factor in eq. (^7]). Thus this term is related to the shape of 
the freeze-out hyper-surface distribution ( which distribution tends to a single hyper-surface 
if At 0). 

The parameters of the BECF-s are dominated by the smaller of the geometrical and the 
thermal scales not only in the spatial directions but in the temporal direction too according 
to eqs. ( ^ - |58|) . These analytic expressions show that even a complete measurement of the 
parameters of the BECF as a function of the mean momentum K may not he sufficient 



to determine uniquely the underlying phase-space distribution [p6| ,p|,p!2|,|9|,p!3| . We also can 
see that the LCMS frame approximately coincides with the LSPS frame for pairs with 
I y^ — Y \« 1 + Arj^Mt/T^ — Arj^ and the terms arising from the non- vanishing values of rjg 
can be neglected. In this approximation, the cross-term generating hyperbolic mixing angle 
Tjs thus we find the leading order LCMS result: 

C{Ak- K) = l + K expi-RlQl - Rl,,Ql,, - i?LQL), (59) 

with a vanishing out-long cross-term, Rout,L = 0. To leading order, the parameters of the 
correlation function are given by 
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Rlide - (60) 

RI, = RI + P^AtI (61) 
Rl = r^^vl (62) 

Observe that the difference of the side and the out radius parameters is dominated by the 
hfetime- |?arameier Ar^,. Thus a vanishing difference between the R^^^ and R'^i^e can be 
generated dynamically if the duration of the particle emission is large, but the thermal 
duration Att becomes sufficiently small, c.f. eq. (|55|). This in turn can be associated with 
intensive changes in the local temperature distribution during the course of the particle 
emission. 

Observe, that the BECF in an arbitrary frame can be obtained from combining eqs. 



581) with the general expressions given by eqs. (p7| - p^ ). In that case, the value of r], 



Y + j^LCMS ^Y + {yo- Y)/{1 + Ar]\l/Ar]^ - 1)) has to be used in eqs. d^H). 



Note also that in our results higher order terms arising from the non-vanishing value of 



rjs in the LCMS are summed up, while in refs. [|T2[ the first sub-leading corrections were 
found. 



VII. INVARIANT MOMENTUM DISTRIBUTIONS 



The IMD plays a complementary role to the measured Bose-Einstein correlation func- 
tion P, p!3| , |26[] . Thus a simultaneous analysis of the Bose-Einstein correlation functions and 
the IMD may reveal information both on the temperature and flow profiles and on the 
geometrical sizes. 

For the considered model, the invariant momentum distribution can be calculated as 
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where the geometrical contribution to the effective temperature is given by TG{mt) = 
ToRl/Rlimt) = (a2 + b'^)MtRl/T^ and the fraction / is defined as / = b^ia^ + b^), 
satisfying < / < 1. In eq. (BSf) the quantities ArjJ^ and r/J are defined as 



1 1 1 

+ cosh[r7l], (64) 



Ts = {yo-y)/{i + AvV^VT), (65) 

i.e. these quantities differ from the Arj and rj^^^^ by the contributions of the Cooper-Fry 
pre-factor. This happens because eq. (|63D is obtained by applying the saddle-point method 
for eq. with the model emission function eq. (0) in such a way that the Cooper-Fry 
pre-factor mtCos\i[ri — y\ is kept exactly and the remaining factors are approximated with 



the saddle-point technique in LCMS described in details in the previous section. Since 
the saddle-point equations are solved in LCMS in a region where rjs—y << 1, a term in the 
exponent is approximated by cosh.[ris — y] ~ 1 + (r^^ — y)^/2. 

For the considered model, the rapidity-width Ay{mt) of the invariant momentum dis- 
tribution at a given rrit, shall be dominated by the longer of the thermal and geometrical 
length-scales. If the condition a ~ is fulfilled, i.e. / ~ 1, the longer of the thermal and 
geometrical scales shall also dominate T*, the effective temperature (slope parameter) at a 
mid-rapidity ?/o- E.g. the following relations hold: 

1 / 1-/ 
Ay^(mt) = Ar/^ + Arj'^^mt) and — = 1 . (66) 

To + TG{mt = m) Tq 

That is why the IMD measurements can be considered to be complementary to the BECF 
data. 

In the special limiting case when gradients of the temperature are negligible, a = 
and / = 1, we have = Tq + mb'^R'Q/TQ . If the flow velocity at the geometrical radius 
(ut) = bRc/TQ is independent of the particle type, we obtain a relation 

T,=TQ + m{utf. (67) 
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A similar relation for the mass dependence of the mean (transverse) kinetic energy was 
introduced in ref. for longitudinally expanding systems with non-relativistic transverse 



flows. This simple relation ( ^) for the effective temperature can be considered as a special 
case of the more general eq. (|66D. 

The measured IMD can be obtained from the IMD of the core as given above and from 
the measured A,f(p) parameters of the BECF as 

1 

Niip) = ^=iVi,c(p), (68) 
'Kip) 



see e.g. ref. ( |^B[) for further details. 

The invariant momentum distribution described by eq. ( |53D features two types of low 
transverse momentum enhancement as compared to a static thermal source with a slope pa- 
rameter T^,. One may introduce the volume factor or K(?/, rrit) which yields the momentum- 
dependent size of the region, where the particles with a given momentum are emitted from: 

K(2/,m,) = (27rA^/r2)i/2(2vri?2)^. (69) 

At 

The rapidity-independent low-pt enhancement is a consequence of the transverse mass de- 
pendence of this effective volume, which may depend on rrit for certain limiting cases in the 
following ways 

fe/2 



(70) 



where k = for a static flreball [a = b = d = and At]^ >> Arf). The case k = 1 
is satisfled for a = b = d = and Ar/^ >> Arj^, which describes long, one dimensionally 
expanding finite systems The case k = 2 corresponds toa = b = Oy^d and Arj'^ >> Arj^, 
describing longitudinally expanding systems with cooling. The case k = 3 corresponds to 
a^O,b^O = d i.e. three-dimensionally expanding, cylindrically symmetric, finite systems 
possibly with a transverse temperature profile p6[ and the k = 4 case corresponds to the 
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same but with a d ^ parameter, describing the temporal changes in the local temperature 
during the particle emission process appended with the condition Att << Ar. Thus the 
inclusion of this effective rrit dependent volume factor into the data analysis not only would 
undoubtedly increase the precision of the measurements of the slope parameters, but in turn 
it also could shed light on the dynamics of the particle emission from such complex systems. 

The rapidity- dependent low-pt enhancement, which is a generic property of the longitudi- 
nally expanding finite systems [^, reveals itself in the rapidity-dependence of the effective 
temperature, defined as the slope of the exponential factors in the IMD in the low-pt limit 
at a given value of the rapidity. The leading order ||29| result is 



^ ^ -2 

^e//(2/) = with a = 



Ar/2 + _ 



(71) 



y m J 



l + a{y- yaf Ivr? 

Please note that this analysis of the low transverse mass region of the IMD relies on the 
applicability of the saddle-point method in the low transverse momentum region too. Thus 
it may be valid for kaons or heavier particles (as well as locally very cold pionic systems). 
However, in case of pions, the self-consistency of the applied formulas and their region of 
validity has to be very carefully checked. This region of the applicability of the saddle-point 
technique for the considered finite systems is discussed in detail in the next section. Note 
that the low transverse momentum region is populated by a number of resonance decays. 



For the long-lived resonances, thus a non-trivial 1/ y A*(p) factor may appear and contribute 
to both the rapidity-dependent and the rapidity- independent low-pt enhancement. Although 
this factor is measurable from the shape- analysis of the BECF, care is required to study the 
contribution of the decay products of short lived resonances to the momentum distribution 
of pions. Kaons or other heavier particles thus provide a cleaner test for these analytic 
results as compared to pions. 

The high-pt enhancement or decrease refers to the change of the effective temperature 
at mid-rapidity with increasing rrit . The large transverse mass limit Too shall be in general 
different from the effective temperature at low pt given by since 
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Too = and — 



( Tcim) \ 

1-/ 

^ To + TG{m)j 



(72) 



2 - / T, 2 - / 

Utilizing Tq/Tq = Rq/B^, the high-pt enhancement or decrease turns out to be controlled 
by the ratio of the thermal radius RT{mt = m) to the geometrical radius Rq- One obtains 
Too > ^* if R^im) > Rq and similarly Too < ^* if Rj'i'^) < Rg- Since for large colliding 
nuclei Rq is expected to increase, a possible high-pt decrease in these reactions may become 
a geometrical effect, a consequence of the large size. 

VIII. LIMITATIONS 

The simple analytic formulas presented in the previous sections are obtained in a saddle- 
point approximation for the evaluation of the space-time integrals. This approximation 
is known to converge to the exact result in the limit the integrated function develops a 
sufficiently narrow peak, i. e. both Arjl << 1 and At^/t^ << 1 are required. This in 
turn gives a lower limit in rrit for the applicability of the formulas for the class of models 
presented in the previous section. 

[However, the saddle-point method may give precise results even when the integrand does 



not develop a narrow peak. For the emission function of ref. |J3[ the saddle-point method 
gives exact result, because that emission function can be re- written in Gaussian form.] 
From the requirement Arj^ << 1 we have 

mtCosh[r7,]/To » 1 + l/cosh2[r7,] - 1/Ar]^. (73) 

From the requirement At^/t^ << 1 one gets 

mt cosh[r/,]/To » (1 - T^/AT^)/d^ (74) 

These are the conditions governing the validity of the calculation of the parameters of the 
BECF as presented in section |VI[ Compared to the condition (|73D, the condition of validity 
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of the calculation of the invariant momentum distribution is less stringent, since one needs 
to satisfy only Aril << 1) which yields 

rrit cosh[rjJ/To » 1 - l/Ar/^. (75) 

For the case of the NA44 measurement one may estimate the region of reliability for the 
analytic formulas presented in the previous parts using inequalities eq. (|73| , |7^ . Since the 
data indicate Rside — Rout within errors, the inequality Ar^/r| << 1 and its consequence, 
eq. (|7^ seems to be well justified. In the inequality eq. ( [73|) the finite longitudinal size 



plays an important role. For infinite systems, Ar/ = oo, the calculations of BECF parameters 
are reliable for nit » 2To (since ?7s = for infinite systems in LCMS), while the calculations 
for the IMD are reliable to >> Tq. In the mid-rapidity region where NA44 data were 
taken, one has ?7s — and for finite systems one finds >> To(2 — 1/ Arf). Note that this 
estimated lower limit in irit is extremely sensitive to the precise value of At] in the region 
Ar/ ~ 1/v^ ^ 0.7. Thus for finite systems the region of applicability of our results extends 
to lower values of uit than for infinite systems which were recently studied in great detail 
in ref. I^O] . The inequality eq. (^) can be used in basically two ways: either one assumes a 



value for Tq and then obtains the lower limit in for the applicability of the saddle-point 
method, or one assumes that the saddle-point method is applicable to a certain value of mt 
(e.g. in case it gives a good description of data), and then one obtains an upper limit for the 
corresponding Tq parameter, the central temperature at the mean time of last interactions. 

An upper transverse momentum limit is obtained for the validity of the calculations in 
sections [V^[VII| from the requirement r^; ,,/to < 1 or tI^/tq << 1. This condition and the 



requirement brx^s/To < 1 are fulfilled simultaneously if 

A « . (76) 

6max(l, b) max(l, a, b) 

If + 6^ ~ 1 this condition simplifies to (3t << 

When comparing to data, detailed numerical studies may be necessary to check the 

26 



precision of the saddle-point integration. In a subsequent paper we plan to present these 
studies together with a detailed comparison of the model to the available NA44 data. 

IX. LIMITING CASES 

Observe that both the thermal and the geometrical length-scales enter both the pa- 
rameters of the Bose-Einstein correlation function and those of the invariant momentum 
distribution. We limit the discussion in this part to the mid-rapidity region rj^^'^'^^ ~ and 
we assume Ai]"^ << 1, i.e. we neglect the 1/ cosh^[r75] term in eq. (p^ . Various limiting 
cases can be obtained as combinations of basically the relative size of the thermal and the 
geometrical scales in the transverse, longitudinal and temporal directions. These in turn 
are: 

i) If RxiMt) >> Rg in a certain Mj interval, we have also Tg >> Tcimt) at the same 
transverse mass scale. In this region, the side radius parameter shall be determined by the 
geometrical size Rside = R* — Rgi hence it shall be transverse mass independent. 

The rrtt distribution at mid-rapidity shall be proportional to exp(— m^/To). 

a) If ArjT » Ai], we have Rl — tqAt] and the rapidity-width of the IMD shall be 
dominated by the thermal scale, Ay^{mt) ~ Arj^ = TQ/rrit. 

in) If Att » At, the temporal duration shall be measured by Rl^^ — Rli^e — Pt^^^- 
The invariant momentum distribution shall be influenced only through the At^/At ^ 1 
factor in V^,. 

These cases are rather conventional limiting cases. An unconventional limit complements 
each: 

iv) If RxiMt) << Rg in a certain Mt interval, we have also Tq << Tcimt) at the 
same transverse mass scale. In this region, the side radius parameter shall be determined 
by the thermal size Rside = -R* — RriMt), hence it shall be transverse mass dependent, 
i?Le « r^To/M,. 

The mt distribution at mid-rapidity shall be proportional to exp(— m^/T^). If << 6^, 
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we have ~ TG{m) as follows from eq. (^). 

v) If ArjT « Ar], we have the leading order LCMS result Rl — tqAtit ~ tq^Tq/Mi^ 
and the rapidity- width of the IMD shall be dominated by the geometrical scale, Ay^ ~ Arf. 

vi) If Att « At, the thermal duration shall be measured by -R^ut ~ R-lide — Pt^Tr — 
P^t^Tq/ {(PMt). For large values of the transverse mass, the model thus shall feature a 
dynamically generated vanishing duration parameter, which has a specific transverse mass 
dependence. The invariant momentum distribution shall be influenced only through the 
At^/At ~ \j ^Jnit factor in K. 

Some combinations of cases i) - vi), are especially interesting, as: 

vii) If all the finite geometrical source sizes, RcAt] and Ar are large compared to the 
corresponding thermal length-scales we have in LCMS 

Ar^ ^ , (77) 

Mt 

Rl ^ — , (78) 



RLe ^ • (79) 



+ Mt 



Thus if (i^ >> + ^ 1 the model may feature a dynamically generated vanishing duration 
parameter. In this case, the model predicts an Mt - scaling for the duration parameter as 



1 

Ar^ (X — . (80) 



Mt 



This prediction could be checked experimentally if the error bars of the measured radius 
parameters were decreased to such a level that the difference between the out and the side 
radius parameters would be significant. 

Alternatively, if the vanishing duration parameter of the BECF is generated due to a 



very fast hadronization process as discussed in ref. then one has 
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At^ ~ Ar^ oc const, (81) 

i.e. in this case the duration parameter becomes independent of the transverse mass. 

If the finite source sizes are large compared to the thermal length-scales and if we also 
have + 6^ ~ 1, one obtains an Mt -scaling for the parameters of the BECF, 

To 1 
Rlde-Rlut-Rl-rl — , valid for A « -. (82) 

Mt b 

Note that this relation is independent of the particle type and has been observed in the 
recent NA44 data in S + Pb reactions at CERN SPS [§. Preliminary NA49 data for Pb + Pb 
at CERN SPS are also compatible with this scaling law [^. This M^-scaling may be valid 
to arbitrarily large transverse masses with ^ 1 if b << 1. The lower limit of the validity 
of this relation is given by the applicability of the saddle-point method, eqs. (^,0). To 
generate a vanishing difference between the side and out radius and an Mt-scaling for the 
BECF radii simultaneously, the parameters have to satisfy the inequalities ( [73|J7^ ) as well 
as b << + 6^ ^ 1 << d"^, i.e. the cooling should be the fastest process, the next dominant 
process within this phenomenological picture has to be the development of the transverse 
temperature profile and finally the transverse flow shall be relatively weak. If the temporal 
changes of the temperature are not intensive enough than a small life-time parameter can 



also be obtained by a fast hadronization and simultaneous freeze-out as discussed in ref. |27 
with At ^ 0. 

We would like to emphasize that there are a number of conditions in the model which 
need to be satisfied simultaneously to get the scaling behavior, which is supported by 9 
NA44 data-points (3 for kaons and 6 for pions, ref. 0). One has to wait for future data 
points to learn more about the experimental status of the scaling. The model presented in 
this paper may describe more complex transverse momentum dependences of the parameters 
of the Bose-Einstein correlation function, too, the Mj-scaling is only one of its virtues in a 
specific limiting case. However, it is rather difficult to get a limiting case with -R^ ~ Rside ~ 
Rout oc 1/ y/Mf in analytically solvable models. Such a behavior is related to the cylindrical 
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symmetry of the emission function. 

Thus the symmetry of the BECF in LCMS can be considered as a strong indication for 
a three- dimensionally expanding, cylindrically symmetric source, possibly with a transverse 
and temporal temperature profile. The LCMS frame is selected if the mean emission point or 
saddle-point stays close to the symmetry axis even for particles with a large transverse mass, 
r^^s(^t) << To and if the finite longitudinal size introduces only small difference between 
the LCMS and LSPS frames, i.e. \y-yo\« 1 + Ar/2(mj/To - 1) . In the case considered in 



section ^ the emission function is cylindrically symmetric and so the BECF is symmetric 



in the LCMS of the pair ( and not in the center of mass system of the pair ||32|| ). 



via) It is interesting to investigate the other limiting case when Rt >> Rg, ^Vt » 
and Att >> At by combining the limiting cases i, - Hi. In this case one obtains 

i?i^roW = <G, Rlde^Rh Rl^^Rl + P^Ar', (83) 

T 

Ay^imt) ~ Ar,l = — , = To. (84) 

Thus, if the thermal length scales are larger than the geometrical sizes in all directions, the 
BECF measurement determines the geometrical sizes properly, and the pt and the dn/dy 
distributions are determined by the temperature of the source. In this case the momentum 
distribution reads as 

( mt cosh.{y - yo)\ 



Ni{p) oc mt cosh(?/ - yo) exp 



(85) 



V ^0 J 

which is a thermal distribution for a static source located at the mid-rapidity yo- 

Thus two length-scales are present in all the three principal directions of three- 
dimensionally expanding systems. The BECF radius parameters are dominated by the 
shorter of the thermal and geometrical length scales. However, the rapidity-width of the 
d^N/dy/dmf distribution, Ay'^init) is the quadratic sum of the geometrical and the thermal 
length scales, thus it is dominated by the longer of the two. Similarly, the effective tem- 
perature is dominated by the higher of the two temperature scales for / ~ 1 according to 
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eq. (|66D. The effective temperature of the mt distribution is decreasing in the target and 
projectile rapidity region in this class of analytically solvable models. 



This study is a generalization of the basic ideas presented and illustrated in ref. |T3 
for the case of three- dimensionally expanding, cylindrically expanding finite systems with a 
scaling longitudinal flow, weak transverse flow and a transverse and temporal temperature 
profile. 



X. SUMMARY 

A general formulation is presented for the two-particle Bose-Einstein correlation function 
for cylindrically symmetric systems undergoing collective hydrodynamic expansion, c.f. eqs. 

. Note that these relations were shown to be valid for certain limited classes 
of emission functions. The resulting class of Bose-Einstein correlation functions, however, 
includes non-Gaussian correlation functions too. The case of Gaussian correlation functions 
is studied in detail and the radius parameters are expressed in the LAB, LCMS and LSPS 
systems, where the functional form of the correlation functions becomes more and more 
simplified. The cross-term generating hyperbolic mixing angle is identified with the value of 
the rj variable of the saddle point in the considered frame. 

A class of Gaussian models is introduced which in some regions of the model-parameters 
may obey an MfScaling for the side, out and longitudinal radius parameters. Vanishing 
effective duration of the particle emission may be generated by the temporal changes of the 
local temperature during the evaporation. The model predicts an Mt-scaling also for the 
duration parameter in this limiting case. 

Finally we stress that both the invariant momentum distribution and the Bose-Einstein 
correlation function may carry only partial information about the phase-space distribution 
of particle emission. However, their simultaneous analysis sheds more light on the dynamics 
and the geometrical source-sizes. 
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APPENDIX 

In this appendix we give a simple example when the Fourier-transformed emission function 
exists but the Gaussian version of the saddle-point method is not applicable. Let us consider 
the one-dimensional Lorentzian distribution function 

1 1 

S{r) = . (86) 

7iR (1 + rVi?') 

Here r is a real variable (in one dimension). The corresponding correlation function is 

C{q) = 1+ I S{q) \\ (87) 



with 



which yields 



/oo 
dr S{r) exp{—iqr), (88) 
-oo 



C(g) = 1 + exp(-2 I g I i?). (89) 



This function is not analytic at q = because it depends on the modulus of q, and for 
positive values of q its Taylor expansion starts with a linear term. This is to be contrasted 
with the results for the saddle-point method. If the Gaussian version of the saddle-point 
method is applicable, then S{q) can be expanded into a Taylor series around g = as 
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S{q) = l + i{r)q- {r^)q^/2 + ... . (90) 
Here the average of a function of variable r is defined as 

/oo 
drf{r)S{r) (91) 
-oo 

and the two-particle correlation function can be written as 

C{q) ^ 2 - q^Rl ^ 1 + exp{-q'^Rl) (92) 

with 

Rl = (r^) - {rf. (93) 

Since for the considered function Rq = oo, the Gaussian saddle-point method is not ap- 
plicable. Still, the Fourier-transformed emission function and the BECF exist, as given by 
eq. (ii. 

Cases similar to this are characterized by non-Gaussian correlation functions Similar 
examples can be found among multi-variate distributions, too. 
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